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Abstract. In this article, we consider a stochastic PDE of parabohc type, driven by a space-time 
white-noise, and its numerical discretization in time with a semi-implicit Euler scheme. When the 
nonlinearity is assumed to be bounded, then a dissipativity assumption is satisfied, which ensures 
that the SDPE admits a unique invariant probability measure, which is ergodic and strongly mixing 
- with exponential convergence to equilibrium. Considering test functions of class , bounded and 
with bounded derivatives, we prove that we can approximate this invariant measure using the 
numerical scheme, with order 1/2 with respect to the time step. 

1. Introduction 

In this paper, we are interested in the discretization in time of the following stochastic reaction- 
diffusion equation 

(1) = + y{t. i)) + 

for t > 0,.^ G (0, 1), with the initial condition y(0,^) = and homogeneous Diriclilet boundary 

conditions y{t,0) = y{t, 1) = 0. The stochastic perturbation ^^g^^^ is a space-time white noise: the 
rigorous interpretation of ([T]) is given by an abstract evolution equation ([2]) - in the sense of [3j - in 
the Hilbert space H = L^(0, 1), driven by a Wiener cylindrical process - see Section [2.31 

The numerical approximation of stochastic equations has been extensively studied during the last 
thirty years. First we recall that one can look at strong approximation results - when the trajectories 
of the continuous and the discrete-time processes are compared - or at weak approximation results 
- when the laws at a fixed time are compared. The simplest method in the case of SDEs is the 
Euler-Maruyama scheme; it is built as a straightforward extension of the well-known explicit Euler 
method for ODEs, which is of order 1: if we consider in a SDE with regular coefficients 

dXt = f{Xt)dt + a{Xt)dBt, Xo = X, 
its numerical approximation is defined for a given step-size At by 

Xq = X, 

Xn+l =Xn + fiXn)At + - BtJ. 

Due to the regularity properties of the Brownian Motion, this scheme is in general only of order 
1/2 in the strong sense - i.e. for nAt < T we have E\XnAt - X„P < C(r)At - while it is of order 
1 in the weak sense - i.e. for test functions (p of class C^, with bounded derivatives, we have a 
bound \E,(p{XnAt) — < C{T,(j))At. The idea for proving that weak order is 1 and not 1/2 

is to consider the Kolmogorov equation associated with the process Xt, which is satified by the 
function {t,x) i— )■ E,(j){X{t,x)) - see [26j, [28j. The books [14j and [22] - see also [21j - contain various 
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numerical schemes - like the well-known Milstein scheme, some implicit schemes, and methods based 
on stochastic Taylor expansions - with their order of convergence in both strong and weak senses. 

Numerical methods for SPDEs like ([1]) need discretization both in time and in space. For example, 
time discretization leads to explicit or implicit methods, while discretization in space can be done 
with finite difference or finite element methods. Basically, the result for space-time white noise 
driven equations is convergence with strong order 1/2 in space and only 1/4 in time - under some 
Courant-Friedrichs-Lewy conditions when necessary: see [5], |10| . |11| . |12| . |13| . |29| . 

If we look at the abstract formulation of ([1]) in the Hilbert space H, results have also been proved 
for the time discretization using semi-implicit Euler schemes: the strong order of convergence is 
1/4 - see [21] - while the weak order of convergence is 1/2 - see [6]. Moreover in ||8j, the authors 
have studied weak convergence in the case of linear equations when using a finite element method 
in space. We follow here the framework of [6]: we consider the stochastic evolution equation in H 
and we use a time discretization with a semi-implicit Euler scheme, with no discretization in space. 

In this work, we are interested in the behaviour of the weak convergence estimates when the final 
time T goes to infinity: can we replace constants C{T,(p) by a constant C(0) independent from 
T? Passing to the limit, we thus ask the more general following question: can we use a numerical 
scheme to approximate the invariant probability measure of the continuous time process - which is 
assumed to be unique, ergodic and with exponential convergence to equilibrium? The SDE case has 
been studied with various methods: in |27| . the weak error analysis is made by showing that the 
time derivatives of the solution of the Kolmogorov equation are exponentially decreasing in time; 
in |17| . some general conditions are given for the ergodicity of the numerical scheme, thanks to the 
theory of geometric ergodicity of Markov Chains. In |19| . the approximation result is shown thanks 
to the use of a Poisson equation. 

In the case of SPDEs, general ergodicity results have only been obtained recently - see Section 
14.31 - and the problem of approximation of invariant measures by numerical schemes has not been 
studied yet. 

We prove the following result: 

Theorem 1.1. For any < k < 1/2, tq > and for any function <j), there exists a constant 
C > such that for any m > 2, y £ H and < r < tq 



The continuous process Y is defined by Equation ([2]) below, and the numerical approximation 
{Ym{T,y)) with time step r and initial condition y is defined by ([8]). 

We notice that the right-hand side of the previous estimate contains a singularity when m = 1, 
which is due to a lack of regularity of the infinite-dimensional processes - details are given in Sub- 
Section 16.11 below. However, if we look at the error at a fixed time T = mr, if r is small enough 
we just need to change the constant C = C{T); moreover since we are interested in the behaviour 
when m goes to infinity, this term plays no role. 

With the assumptions precised below, we show in Section 14.31 that the SPDE admits a unique 
invariant probability measure /I, which is ergodic and strongly mixing, with exponential convergence 
to equilibrium; nevertheless we can in general only show the existence, not the uniqueness, of 
invariant measures for the numerical approximation, and we can prove the following result: 

Corollary 1.2. For any 0<k<1/2, ro>0 and for any function (f), there exists constants 
c > 0, C > such that for any < r < tq, any initial condition y £ H and any m > 1 



n4>{Y{mT,y))] - n4>{YUr,ym < C{1 + |y|3)(((m - l)r)-i/2+- + l)r 



1/2-K 



n4>{Ym{r,y))]- [ (t>dn<C{l + \y\^){ 



m 



1/2-K 



1 



+ T 



/2-) + C(l + |y|2)e 




Moreover, if /j,'^ is an ergodic invariant probability measure of the numerical scheme (5^(t, .))meN; 
we have 

I / (pdfi^ - [ (pd4l\< C7ri/2-«. 
Jh Jh 

all the ergodic invariant measures of the numerical approximation are then close to the unique 
invariant probability measure of the continuous process. 

Up to our knowledge, this is the first result of this kind for SPDEs. 

The key point, like in |27| . for obtaining bounds independent from the time T = mr, is to prove 
that derivatives of the solution u of the underlying Kolomogorov equation - mentioned above - 
decrease exponentially in time: this is done in Section 15.21 using a coupling method. 

Moreover, an essential tool in [6] is the use of Malliavin calculus and of an integration by parts 
formula in order to transform a stochastic integral - which is not regular enough in space in the 
infinite dimensional setting - into an expression which can be controlled; in Lemma [4. 5| we see that 
the involved Malliavin derivatives may increase exponentially fast with respect to time. The solution 
is to separate the lack of regularity problem and this badly controlled growth with a decomposition 
of the interval into two parts: in the first one we need the integration by parts formula, and we can 
control the Malliavin derivatives, while in the other one we can directly give an appropriate bound. 
We also provide an improvement with respect to [6j: we can consider a more general nonlinear 
coefficient G - like a Nemytskii operator, see Example 12.71 

We also notice that some more general equations, with additive noise which is white in time but 
colored in space, can be studied with our method, with suitable assumptions on the coefficients. 
For a very smooth noise, the numerical analysis on finite time is easier to treat, but then ergodic 
properties and long-time behaviour require different techniques; this will be treated elsewhere. 

The case of some equations with multiplicative noise which satisfy the Strong Feller Property is 
also covered by our technique of proof, but with some additional difficulties - see for instance the 
restrictive condition on the diffusion coefficient in [6J. Since all the necessary ideas are contained 
here and in [^, we only focus on the additive noise case. 

The paper is organized as follows: in Section [21 we precise the assumptions made on the coefficients 
of the equations, and we define the numerical method in Section [3l In Section IU we give some 
bounds on the solutions of the continuous and discrete equations, and we study their asymptotic 
behaviour: existence of invariant measures, and uniqueness for the continuous equation, following 
from Proposition 14.71 In Section [5l we explain the proof of Theorem II. H and we give a proof of 
Corollary II. 2j more precisely, in Section 15.21 we prove the exponential decreasing in time of the 
derivatives of the function u. Eventually Section [6] contains the proofs of the remaining estimates 

2. Notations and assumptions 

Let H he a, separable Hilbert space, with norm denoted by \.\h or simply |.|. We consider 
equations of the form 

dY{t, y) = {BY{t, y) + G{Y{t, y)))dt + dW{t) 

In the next paragraphs, we explain the assumptions on the linear operator B and on the nonlinear 
coefficient G; we also recall how the cylindrical Wiener process W is defined, and how we can 
construct solutions of this equation. 

2.1. Test functions. To quantify the weak approximation, we use test functions (p in the space 
C^(F,M) of functions from F to M that are twice continuously differentiable, bounded, with first 
and second order bounded derivatives. 
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Remark 2.1. In the sequel, we often identify the first derivative D(j)(x) £ C{H,M) with the gradient 
in the Hilbert space H, and the second derivative D'^(j){x) with a linear operator on H , via the 
formulas: 

< D(j){x)^ h >= D(j){x).h for every h £ H 
< D^(l){x).h, k >= k) for every h,k £ H. 

In the sequel, we use the following notations: 

||$||oo = sup \^{x)\h 

W^Wi = sup \D^{x)\h 
||$||2 = sup \D^<^{x)\c{H)- 

2.2. Assumptions on the coefficients. 

2.2.1. The linear operator. We denote by N = {0, 1,2,.. .} the set of nonnegative integers. 
We suppose that the following properties are satisfied: 

Assumptions 2.2. (1) We assume that there exists a complete orthonormal system of elements 
of H denoted by {fk)ken, md a non- decreasing sequence of real positive numbers {nk)k£N 
such that: 

Bfk = -l^kfk for all k eN. 
(2) The sequence (fik) goes to +oo and 

+00 ^ 

E— < +00 ^ a> 1/2. 
11°' 

k=0 

The smallest eigenvalue of — i? is then hq. 

Example 2.3. The equation ^ enters in this framework: we can choose B = with the domain 
H'^{0,1) n//Q(0, 1) C L^(0,1) - corresponding to homogeneous Dirichlet boundary conditions. In 
this case for any k £ 'N Hk = 7r^(/c + 1)^, and fk{^) = V^sm{{k + l)7r^) - see [2\. 

For a G {1, 2, . . .}, we denote by H]sf the subspace of H spanned by /o, . . . , /tv-i, and by P/v 
the orthogonal projection of H onto H^. 

The domain 0(3) of B is equal to D{B) = {y = Yl=oykfk S H,Yl=o{lJ^k?\yk? < +oo}. We 
can more generally define fractional powers of —B, for b G [0, 1]: 



{-Bfy = Y,^jiykfk(^H, 

k=0 

with the domains 

{+00 +00 
y = Yl ykfk e H, \y\l := Y.il'kfV < +00 
fc=0 k=0 



When b G [0,1], we can also define the spaces D{—B) and operators {—B) , with norm 

denoted by when y = J2k=oykfk G H, we have {-B)-''y = YJk=oH^ykfk and \yW := 

^-2^)u, |2 
2^k=oif^k) ml ■ 



The semi-group {e^^)t>o can be defined by the Hille-Yosida Theorem - see [2]. We use the 
following spectral formula: if y = Ylt^ Ukfk ^ H, then for any t > 

+00 

k=0 

For any t > 0, e*'^ is a continuous linear operator in H, with operator norm e~^°^. The semi-group 
(e*'^) is used to define the solution Z{t) = e^^ z of the linear Cauchy problem 

'^^ = BZ{t) with Z(0) = z. 

To define solutions of more general PDEs of parabolic type, we use mild formulation, and Duhamel 
principle. 

This semi-group enjoys some smoothing properties that we often use in this work. Basically we 
need the following properties, which are easily proved using the above spectral properties. 

Proposition 2.4. Under Assumption \2.SX for any a G [0, 1], there exists > such that we have: 

(1) for any t > and y & H , 

|e*^2/U < C^t-''e-'^'\y\H. 

(2) for any < s < t and y ^ H , 

W"y-e'^y\H <C„- — -^e 2 ' y 

(3) for any < s <t and y G D{-Ay , 

|e*^y - e^^yl^, < a(t - sfe-'^lyU. 

2.2.2. The nonlinear operator. The nonlinear operator G is assumed to satisfy some general assump- 
tions. In Example 12. 7| we give the two main kind of operators that can be used in our framework. 

Assumptions 2.5. The function G : H ^ H is assumed to he hounded and Lipschitz continuous. 
We denote hy Lq the Lipschitz constant of G 

We also define for each N > 1 a function Gn ■ Hn Hn, with GN{y) = PNG{y) for any 
y G Hn- We assume that each Gjy is twice differentiahle, and that we have the following hounds on 
the derivatives, uniformly with respect to N : 

• There exists a constant Ci such that for any N > 1, y £ and h G Hn 

\DGN{y)MH<Ci\h\H. 

• There exists r] G [0, 1) and a constant G2 such that for any N > 1, y £ Hn and any 
h,k £ Hn we have 

\{-BrD^GN{y).{h,k)\ < C2\h\H\k\H. 

• Moreover, there exists a constant C3 such that for any N > 1, y £ Hn and any h,k £ Hn 

\D^GN{y).ih,k)\ <C3\h\^_B)v\k\H. 

Since G is bounded, the following property is easily satisfied: 

Proposition 2.6 (Dissipativity). There exist c > and C > such that for any y £ D{B) 

(3) <By + G{y),y><-c\y\^ + C. 

We remark that we have uniform control with respect to the dimension of the bounds on Gtv 
and on its derivatives, and that ([3]) is also satisfied for Gat, with constants c and C independent 
from A^. 
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Example 2.7. We give some fundamental examples of nonlinearities for which the previous as- 
sumptions are satisfied: 

• A function G : H ^ H of class C^, bounded and with bounded derivatives, fits in the 
framework, with the choice r/ = 0. 

• The function G can be a Nemytskii operator: let g : (0, 1) x R — M 6e a measurable, 
bounded, function such that for almost every ^ G (0, 1) .) is twice continuously differen- 
tiable, with uniformly bounded derivatives. Then G{y) is defined for every y £ H = L'^{0, 1) 
by 

G{ym=9{^M0)- 

In general, such functions are not Frechet differentiable, but only Gateaux differentiable, with 
the following expressions: 

[DG{y).hm = ^{iMm{o 
[D'G{y).{h,k)m = ^{cMO)Hom- 

If h and k are only I? functions, D^G{y).{h, k) may only be ; however if h or k is L°° , 
it is I?. The conditions in Assumption \2. 5\ are then satisfied as soon as there exists rj < 1 
such that D{—B)'^ is continuously embedded into L°°(0, 1) - it is the case for B given in 
Example \2.'A with r] > 1/4. Then the finite dimensional spaces are subspaces of L°° , 
and differentiability can be shown. 

2.3. The cylindrical Wiener process and stochastic integration in H. In this section, we 
recall the definition of the cyhndrical Wiener process and of stochastic integral on a separable 
Hilbert space H with norm \.\h- For more details, see |3|. 

We first fix a filtered probability space {^,J-, (J-i)t>o, IP). A cylindrical Wiener process on H is 
defined with two elements: 

• a complete orthonormal system of H, denoted by {qi)i^j, where / is a subset of N; 

• a family (/3j)jg/ of independent real Wiener processes with respect to the filtration {{J-t)t>o)] 

then W is defined by 
(4) W{t) = Y,mQ^■ 

When J is a finite set, we recover the usual definition of Wiener processes in the finite dimensional 
space RI^L However the subject here is the study of some Stochastic Partial Differential Equations, 
so that in the sequel the underlying Hilbert space H is infinite dimensional; for instance when 
H = L^(0, 1), an example of complete orthonormal system is (g^.) = (-v/2 sin(A;7r.))fc>i - see Example 
[231 

A fundamental remark is that the series in (jj]) does not converge in H; but if a linear operator 

: H ^ K is Hilbert-Schmidt, then "^W{t) converges in L'^{^l,H) for any t >0. 

We recall that a bounded linear operator ^ : H ^ K is said to be Hilbert-Schmidt when 

k=0 

where the definition is independent of the choice of the orthonormal basis (g^) of H. The space of 
Hilbert-Schmidt operators from H to K is denoted C2{H,K); endowed with the norm \-\c2{H,K) it 
is an Hilbert space. 

6 



The stochastic integral Jg* ^'(s)(iW^(s) is defined in K for predictible processes \I' with values in 

C2{H,K) such that l"^ {^)\'c2{h K)^^ ^ ^-^i moreover when ^ G L^(r2 x [0, t]; £2(-f^, -f^^)), the 
following two properties hold: 

^{s)\'c2{H,K)ds (Ito isometry), 



E [ 'i'{s)dW{s) = 0. 
Jo 



I <v,dW{s)>=^l3k{t)vk; 



A generalization of Ito formula also holds - see [3] . 
For instance, if t; = ^^.gj^ VkQk S H, we can define 

ft 

<W{t),v>-- 

ken 

we then have the following space-time white noise property 

E < W{t),Vl >< W{s),V2 >= t A S < Vl,V2 > . 

Therefore to be able to integrate a process with respect to W requires some strong properties 
on the integrand; in our SPDE setting, the Hilbert-Schmidt properties follow from the assumptions 
made on the linear coefficients of the equations. 

Thanks to Assumption I2.2| it is easy to show that the following stochastic integral is well-defined 
in H, for any t > 0: 

(5) W^{t)= f e^'-'^^dW{s). 

Jo 

It is called a stochastic convolution, and it is the unique mild solution of 

dZ{t) = BZ{t)dt + dWit) with Z(0) = 0. 

Under the second condition of Assumption \2.2\ there exists 6 > such that for any t > we 
have js\G^^\'c^(^H'^ds < -\-oo; it can then be proved that has continuous trajectories - via the 
factorization method, see [3] - and that for any 1 < p < +oo 

(6) Esup|l^^(i)|^ < +00. 

t>o 

We can now define solutions to Equation ([2]), thanks to the assumptions made on the coefficients: 
the following result is classical - see [3]: 

Proposition 2.8. For every T > 0, y G H , the equation ([2]) admits a unique mild solution Y G 
L\Q,C{[0,T],H)): 



(7) Y{t) = e 



tBy^ f eii-^)BG{Y{s))ds+ f e^'-'^^dW{s). 
Jo Jo 



3. Definition of the numerical scheme 
We now define the numerical approximation of Y: denoting by r the time step, we have 
Yk+i{T, y) = Yk{T, y) + TBYk+i{T, y) + TG{Yk{T, y)) + V^Xk+i 

yo{T,y) = y, 

where Xk+i = ^{W{{k + l)r) - W{kT)). 

To simplify the equations, we omit the dependence of on the time-step r and on the initial 
condition y. 
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This expression does not make sense in H. Defining R,- = {I — tB) ^, this last equation can be 
replaced by 

(8) = RrYk + TRrG{Yk) + ^/TRrXk+l. 

which is valid, since Rt is a Hilbert-Schmidt operator on H. 
Remark 3.1. Later, we often use the following expression for Yk: 

k-l k-l 

(9) n = + r ^ R';-'G{Yi) + 

;=o 1=0 

The following expression is also useful: 

k-l ,t 



(10) ^Y.^r-'XI+1 = [ " R'^-'^dW{s), 

1=0 

where Is = [|J - with the notation [.J for the integer part. 

We need the following technical estimate: 
Lemma 3.2. For any < k < 1 and j > I, 

\i-Br^Ri\ciH) < -j-^.jrrhw- 

Proof For any z E H, 



+ 00 -. 



i=0 

+ 00 



" (jr)2(i-) ^ Vl + (1 + //ot)2J'^'^'' 

2 1 1 

4. Preliminary results 

We warn the reader that constants may vary from line to line during the proofs, and that in order 
to use lighter notations we usually forget to mention dependence on the parameters. We use the 
generic notation C for such constants. 

We fix the time step r, as well as m G N; we then introduce the notation T = mT. We also define 
tk = kr. «; > is a parameter, which is be supposed to be small enough. We also control r: for 
some To > 0, r < tq. 

4.1. Galerkin approximation. The first step of the proof is to consider finite dimensional approx- 
imations of the i?- valued processes (^(i))teiR+ ^'^d (5^A;)ifceN: if we fix AT > 1, we define iX^^\t))^^-g+ 
and (Y'^^^^)fegN by the equations 

dFW(t) = BY^^\t)dt + GN{Y^^\t))dt + dW^^\t) 

and 



with the initial conditions Y^^q = Y^^^ = PnU- 

The projection P/v and the nonhnear coefficient Gat have been defined above. VF(^) = PnW 
is a A^-dimensional Wiener process on the subspace Hjsf. We remark that the above equations are 
well-defined on Hj^ - which is a stable subspace of B. 

The important and not difficult to prove result is the following: for any fixed t G and /c G N, 
when N — t- +cxd we have 

E|y(t) - ^ and E|Yfc - Y^^^^ 0. 

We need test functions adapted to the finite-dimensional approximation: for any A'^ > 1, by 
restriction we define $Ar(y) = ^{y) for any y G H]\j] we obtain the following decomposition 

E$(y(mr)) - E$(y„) = E^>(y(mr)) - E$(y(^)(mr)) 

+ E$^(yW(mr)) -E^jvlll"^)) 

+ E$(yif))-E$(y„); 

the first and the third terms converge to when — )• -|-oo. In the sequel, we prove an estimate of 
the second term, which is uniform with respect to dimension A^; letting N — )• +oo then yields an 
estimate on the left hand side. 

Hence we work with the finite dimensional approximation, but we omit the parameter A^. The 
constants appearing below are independent of A^. 

In Section 14.21 we prove some estimates on Y[t) and 1^, and in Section 14.31 we focus on the 
asymptotic behaviour of the processes. 

4.2. Some useful estimates. Bounds on moments of 1^ and 1^ can be proved, uniformly with 
respect to time. 

Lemma 4.1. For any p>l, there exists a constant Cp > such that for every t >0 and y €z H 

E|y(t,y)r<Cp(l + |yn. 
Proof If we define Z{t) = Y{t) - W^{t), we have Z{0) = Y{0) = y, and 

^ = BZ{t) + G{Ym 

and by Proposition 12.61 

^^^=< BZ{t) + G{YmZ{t)> 

=< BZ{t) + G{Z{t)), Z{t) > + < G{Y{t)) - G{Z{t)), Z{t) > 

< -c|Z(t)|2 + C7+ ||G||oo|^(t)| 

< -c'|Z(t)|2+C', 

for some new constants c', G' . 

Then almost surely we have for any t >0 

\Z{t)\<G{l + \y\). 

Thanks to ([6]), the conclusion easily follows. □ 

Lemma 4.2. For any p > 1, tq > 0, there exists a constant C > such that for every < r < tq, 
k and y £ H 

Einr <c(i + i2/n. 
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Proof As in the proof of Lemma [4. II above, we introduce = — Wm, where the process {wm) 
is the numerical approximation of with the numerical scheme ([S]) - with G = 0; it is defined by 

Xm+1 ■ 

Using Theorem 3.2 of [24], giving the strong order 1/4 for the numerical scheme - when the initial 
condition is 0, with no nonlinear coefficient, with a constant diffusion term and under the assump- 
tions made here - we obtain the following estimate: for any p > 1, tq > and < r < 1/2 there 
exists C > such that for any < r < tq and ni > 

(11) Elwrr. - W^(mT)\^P < C7r(i/2-r-)p. 



Thanks to ([6]) and pTj) . we get that for any tq > 0, there exists C > such that for < r < tq and 
m > 

(12) ^Wm\'^ < C. 

Now Zm defined above satisfies Zq = Yq = y and 

Zrn+l = RrZm + 'TRrGiXm)] 

since \Rt\c{H) ^ T+^t^rF' obtain the almost sure estimates 

\Zm+l\ ^ :; — ; \Zm\ + Ct 

1 + A^OT" 

and 

\Zra\<C{l + \y\). 

Thanks to (fT2|) . we therefore obtain the result. □ 
We now introduce the following process: for < A; < m — 1 and tk < t < tk+i 

(13) Y{t) = Yk+ I [BrYk + RrG{Yk)]ds+ [ RrdW{s), 

Jtk Jtk 

where B-t = BRr- The process y is a natural interpolation of the numerical solution (Y^) defined 
by m- Y{tk) = Yk. 

Thanks to Lemma [4. 2 ( we get 

Lemma 4.3. For any p > 1, tq > 0, there exists C > such that for any < r < tq, t > and 
yGH 

E\Y{t)\P < C{l + \y\P). 

In the next Lemma, we give a control on Malliavin derivatives of used in the proof. For an 
introduction to Malliavin calculus, see [23], |25|. The notations here are the same as in [6j, where 
the following useful integration by parts formula is given - see Lemma 2.1 therein: 

Lemma 4.4. For any F G B^''^{H), u € Cl{H) and ^ £ L'^{VL x [0, T], £2(-f^)) an adapted process, 

(14) E[Du{F). ^{s)dW{s)]=E[ Tr(^(s)*L>\(F)D,F)ds], 

Jo Jo 

where DgF : h £ H D^F £ H stands for the Malliavin derivative of F, and B-'^'^(ff) is the set of 
H -valued random variables F = Yli^N ^ifi' '"'^^^ ^ ^^'^ domain of the Malliavin derivative 
for W'Valued random variables for any i, and Yliem Jq '^\DsFi\'^ds < +oo. 

Without any further dissipativity assumption, we are not able to give a uniform control with 
respect to time of the Malliavin derivative of Y. In the proof below, we take care of this problem 
by using this derivatives only at times tk = kr and s such that tk-i^ < 1. 
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Lemma 4.5. For any < /3 < 1 and tq > 0, there exists a constant C > such that for every 
heH,k>l,0<T<To and s G [0, tk] 

\D':Yk\p<C{l + LGT)''-''{l+ ^ 



(l + /.or)(i-«('=-'^)tf 



k la 



Proof For any k > 1, h £ H and s £ [0,tk], using the chain rule for Mahiavin calculus and 
expressions ([9]) and (fTOl) . we have 

k-l 

D^Yk = R'^-^'h + T Rr'^DG{Yi).D'^Yi. 

Indeed, recall that Ig denotes the integer part of -, so that when I < Ig we have D^Yi = 0. 
As a consequence, for k > Ig + 1 

\D'^Yk\ < {l + LGT)''~^'\h\. 

Now using Lemma [3.2| we have 

li-BfD'^Ykl < \ \h\ 



To conclude, we see that 

< C < +00, 

when < T < tq. □ 

4.3. Asymptotic behaviour of the processes. The results of this section are obtained for the 
initial H-valued processes, and for their finite dimensional approximations. 

First, we focus on the existence of invariant measures for the continuous and discrete time pro- 
cesses. We use the well-known Krylov-Bogoliubov criterion - see [4J. Tightness comes from two facts: 
D{—B)^ is compactly embedded in H when 7 > 0, and when 7 < 1/4 we can control moments: 

Lemma 4.6. For any < 7 < 1/4, r > and any y £ H, there exists C{'y,T,y),C{'y,y) > such 
that for any m > 1 and t > 1 

E|y„(r, < C(7, r, y) and K\Y{t, < ^(7, y) 

Uniqueness of the invariant probability measure for the continuous time process (5^(i))teM+ can be 
deduced from the well-known Doob Theorem - see Indeed, since in equation ([2|) noise is additive 
and non-degenerate, the Strong Feller property - see also Lemma 15.41 below - and irreducibility can 
be easily proved. In the proof of the main Theorem II. H we also need speed of convergence, and 
thanks to a coupling argument we get the following exponential convergence result - for a proof see 
Section 6.1 in [7]: 

Proposition 4.7. There exist c > 0, C > such that for any bounded test function (f), any t > 
and any 2/1,1/2 G H 

(15) |E0(y(t,yi)) - E0(y(t,2/2))| < Cm^l + jyip + \y2\^)e-''\ 
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The idea of coupling relies on the following formula: if z/i and U2 are two probability measures 
on a state space S, their total variation distance satisfies 

drv(z^i,i^2) = inf{P(Xi 

which is an infimum over random variables {Xi^X2) defined on a same probability space, and such 
that Xi ~ vi and X2 ~ f2- 

Roughly speaking, the principle is to define a coupling yi, 2/2)1 ^2(^5 2/1, y2))t>o for the pro- 

cesses iY[t, yi)t>o cind Y{[t, y2))t>o such that the coupling time T of Zi and Z2 - i.e. the first time 
the processes are equal - has an exponentially decreasing tail. 

This technique was first used in the study of the asymptotic behaviour of Markov chains - see 
[2]; |16| . |20) - and was later adapted for SDEs and more recently for SPDEs - see for instance 

US], m- 

In fact, uniqueness of an invariant probability measure is an easy consequence of this Proposi- 
tion, and moreover we get for any y £ H and any t >0 

(16) \E<p{Yit,y))- [ cPdJl\<C\\4>\Ul + \y\'')e-'^'. 

Jh 

In general, we do not know whether uniqueness also holds for the numerical approximation 

Remark 4.8. A sufficient condition for the uniqueness of the invariant probability measure of the 
discrete time process {Yk)k£f>s is the strict dissipativity assumption 

Lg < fJ'O, 

where we recall that Lq denotes the Lipschitz constant of G. 

Then trajectories of the processes (lt)tgiK+ and (yfc)fcgN issued from different initial conditions yi 
and y2 and driven by the same noise process are exponentially close when time increases: for any 
tq > 0, there exists c > such that for any < t < tq, k > and t >0 we have almost surely 

\Y{t,yi) - Y{t,y2)\ < e-(^«-^«)*|yi - 2/2! 
\Yk{T,yi) - Yk{T, 2/2)1 < e-^^^|2/i - y2\. 
Proof of uniqueness is then straightforward - and we do not need Proposition ^ - ?| above. 

5. Presentation of the proof of the weak approximation result 

The proof of Theorem II .11 is very technical, so for pedagogy we first introduce the decomposition 
of the error, and identify the term which we control later in Section [6l Some crucial estimates on 
the derivatives of the semi-group with respect to the initial conditions - regularization, long-time 
behaviour - are proved below in Sub-Section 15.21 

5.1. Strategy. We define 

(17) u{t,y)=n<>{Y{t,y))\, 

which is solution of a finite dimensional Kolmogorov equation associated with the finite dimensional 
approximation of ([2]): 

^(i, y) = Lu{t, y) = ^Tr {D\{t, y)) + < By + G{y), Du{t, y) > . 

As explained in the introduction, this is one of the essential tools in the proof of the weak approxi- 
mation result. 
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The weak error at time T = mr can be decomposed with a telescoping sum - where to simphfy 
the dependence of the numerical approximation in r and y is not written: 

E[(^(y(r, y))] - E[(t>{Ym{T, y))] = u{T, y) - E[n(0, Y^{t, y))\ 

m—l 

= V {E[uiT - tk,Yk)] - E[u{T - tk+i,Yk+i)]) 

(18) t^o 

m—l 

= u{T, y) - E[n(r - r, Yi (r, y))] + ^ (a^ + hk + Ck), 

k=l 

where for 1 < k < m — 1 

Ofc = E / < BY{t) - BrYk, Du{T - t, Y{t)) > dt, 
Jtk 

(19) bk=E [ < GiY{t)) - RrG{Yk),Du{T - t,Y{t)) >dt, 

Cfc = -E / Tr((/ - RrR*)D^u{T - t, Y{t)))dt. 

This follows from the use of the Kolmogorov equation and of the Ito formula. 

5.2. Bounds on the derivatives of the transition semi-group. By (fT7|) . u{t, y) = E[0(y(t, y))]; 
since (j) is of class C^, bounded and with bounded derivatives, we are able to prove that with respect 
to y the function u is twice differentiable, and that the derivatives can be calculated in the following 
way: 

• For any /i G H, we have 

(20) Du{t, y).h = E[DcPiYit, y)).r]^^y (t)], 
where ry'*'^ is the solution of 

^a!^ = Br,^^y{t) + DG{Y{t, y)).^^^y{tl 

i]^'y{o) = h. 

• For any h,k £ H, we have 

(21) DMt,y)-ih,k) = E[D''<PiYit,y)).{v''^yit),rj'''y{t)) + D<p{Yit,y)).C''^''^^ 
where (.^'^^y is the solution of 

^C;:^ = i?c"'''^(t) +I)G(y(t,y)).C"''='^(t) + Z)2G(y(t,y)).(ry'^.?^(t),r/'=.?'(t)), 

(^h,k,y^Q^ = 0. 

In [6], the key point for obtaining the weak order 1/2 is to control the derivatives \Du{t,y)\i3 and 
\{-B)'^D'^u{t,y){-B)'y\c{H'j, with (3 < 1/2 and 7 < 1/2 - with the identification of Remark EH 
Moreover, to obtain a long-time weak estimate we need to prove some exponential decreasing of 
such quantities when time t goes to infinity. The two Propositions below are the essential results 
we thus need. 

Proposition 5.1. There exists a constant jl > such that for any < /3 < 1/2, for any t > and 

yeH 

(22) \Duit,y)\p < Cpil + ^)e-^'*(l + |yp). 
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Proposition 5.2. There exists a constant ft > such that for any < /3,7 < 1/2, for any t > 
and y £ H 

(23) \{-BfDMt,y){-By\ciH) < Cp,,{l + i + J^)e-f'\l + Ivl')- 

The singularity t^^ in (j23p is a consequence of tlie regularity properties satisfied by G. Since in 
general during the proof of Tlieorein ll.il we need /3 + 7 to be close to 1, and therefore greater than 
T], only the second singularity t~^~"' plays a role. 

The proofs require several steps. First in Lemma [5. 31 below we prove estimates for a finite horizon 
and general < /3,7 < 1/2; then in Lemma [5.41 we study the long-time behaviour in the particular 
case /3 = 7 = 0; we finally conclude with the proofs of Propositions 15.11 and 15.21 

First, we prove estimates of these quantities for < t < 1 - see Lemmas 4.4 and 4.5 in [6], with 
a difference coming from the assumptions made on the nonlinear coefficient G: 

Lemma 5.3. For any < /3 < 1/2, < 7 < 1/2, there exist constants Gp and C^,7 such that for 
any y £ H and any < t < 1 

\Duit,y)y<^ 

\{-BfDMt,y){-By\ciH) < C^,,{^ + ^). 

If we take another time interval ]0, Tmax] instead of ]0, 1], the constants and C^,^ are a priori 
exponentially increasing in Tmax- 

Proof Owing to (pO|) and (j2ip . we only need to prove the following almost sure estimates, for 
some constants G^ and C^,-y - which may vary from line to line below: for any < t < 1 



(24) 



where the parameter r] is defined in Assumption 12.51 

We use mild formulations, and the regularization properties of the semi- group given in Proposition 



/'?^(i)| = |e*^/i+ I e'^^~'^^DG{Y{s,y)).r]^'y{s)ds\ 
Jo 

<^\h\_p + G j\r^^^y{s)\ds, 



and by the Gronwall Lemma we get the result. 

For the second-order derivative, we moreover use the properties of G in Assumption 12.51 to get 

|^h,fc,y(i)| = I DG{Y{s,y)).C^''''y{s)ds 



+ f e(*-^)^D2^(y(s, y)).{rl'^y{s),r^'''y{s))ds\ 



<G f\C''''^y{s)\ds+ f ^in^\r^^^y(s)\\v''\s)\ds 

Jo Jo [t-sr 

<G j'^ \Q^^''y{s)\ds + Gp,,\h\^p\k\^pt'-^-P-' j'^ ^^—L-^d.. 
To conclude, it remains to use the Gronwall Lemma. □ 
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Thanks to the dissipativity property expressed in Proposition 12.61 we can prove the result in the 
case /3 = 7 = 0. We notice that the proof would be straightforward under a strict dissipativity 
assumption - since then r/'*'^(t) and (^^'^'^{t) would decrease exponentially in t; in the general case 
r]^'y{t) and C^'^'y[t) are exponentially increasing in time so that we can not work directly. Here the 
result comes from the estimate (jlSp of Proposition |47 



Lemma 5.4. There exist constants C and c > such that for any t > and any y £ H 

(25) \Duit,y)\<Ce-'\l + \y\^) and |Z)\(i, < Ce-^*(1 + + 

Proof The Bismut-Elwortliy-Li formula states that if <I> : — )■ M is a function of class with 
bounded derivatives and with at most quadratic growth - i.e. there exists M(<I>) > such that for 
any y £ H we have |$(y)| < M(<^)(1 + - then we can calculate the first and the second order 
derivatives of {t,y) i— )• v{t,y) := 'K^(Y{t,y)) with respect to y. First, we have for any y £ H and 

Dv{t,y).h = f < r,'''y{s),dW{s) > HY{t,y))] 

(26) 

= -E[ / < r^'''y{s),dW{s) > v{t/2,Y{t/2,ym 



t 

the second equality is a consequence of the identity v{t,y) = E,v{t/2,Y{t/2,y)), thanks to the 
Markov property. Using the second formula of (j26p . we obtain a formula for the second order 
derivative: for any y £ H and h,k £ H, 

2 r*/^ 

DMt, y)-{h, k) = -E[ < C'''''yis),dW{s) > v{t/2, Y{T/2, y))] 

(27) 

+ -E[ / < r]^^y{s),dWis) > Dvit/2, y(t/2)).r/'='^(t/2)]. 
t Jo 

We then see, using Lemmas 14.11 and 15.31 - with /? = 7 = - that there exists C > such that for 
any < t < 1, y £ H , h,k £ H 

\Dv{t,y).h\ < ^Mmi + \y\^)\h\, 

(28) 

\DMt,y)-{h,k)\ < jM{<^){l + \y\^)\h\\k\. 

Now when t > I the Markov property implies that u{t, y) = E,u{t — 1, Y{1, y)), and by (|16p we have 

\u{t -l,y)- [ <pdfi\ < Ce-'^'-'\l + |y|2). 

If we choose ^t{y) = u{t -l,y) - (j)dfi, we have u{t, y) = E<I)j(y(l, y)) + ^d/i, with M{(^t) < 
^g_c{t_i)_ ^-^j^ ^ ^-j^g -j^^ obtain for t > 1 

\Du{t,y).h\<Ce'<'-^\l + \y\^)\h\ 

\DMt.y)-{hM<Ce-'^'-^\l + \y\^)\h\\k\. 

Moreover by Lemma 15.31 we have a control when < t < 1 , so that with a change of constants we 
get the result. □ 
We can finally prove the Propositions 15.11 and 15.21 The key tool is the Markov property of the 
process Y which yields the following formula: for any t > 1 

(29) u(t,y) =E[n(t-l,yi(y))]. 
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To get the exponential decreasing, we use Lemma 15.41 at time t — 1 when t > 1, while \h\-p 
appears from rjh^y{l), and with estimates coming from Lemma [5.31 

Proof of Propositions 15.11 and 15.21 Using equation (j29p and Lemma 15. 4| for any i > 1 we have 

\Du{t,y).h\ < Ce'<'-Ma + 1^(1, l')|')|r?'^'^(l)|] < Ce-<'-'\l + \y\')\h\.p, 

where the last estimate comes from Lemmas 14.11 and 15.31 Combining this estimate with the result 
of Lemma 15.31 which gives an estimate for t < 1, we obtain (j22p . For the second order derivatives, 
Lemma 15.31 gives an estimate for t < 1, and for t > 1 we use (|29p to see that 

D\{t,y).{h,k) = E[D^[u{t - l,Y{l,y))].{h,k)] 

= EDMt - 1, ^(1, y)).(r?"'ni), ^?'^'^(1)) + EDu{t - 1, ¥(1, y)).C"''='^(l). 

Using Lemma I5.4| we get an exponential decreasing; thanks to Lemma 14.11 and to the estimates in 
the proof of Lemma 15.31 at time 1 , we obtain 

\DMt,y)-{h,k)\<Ce-<'-'\l + \y\^)\h\.i,\k\.^. 

Then (I23p easily follows. □ 



5.3. Proof of Corollary 11.21 The first estimate is a simple consequence of the Theorem, and of 
the exponential convergence to equilibrium of the continuous-time process - see (fT6|) . We then get 

miYrnir,y))]- [ 0d7l|<C(l + |y|3)(^^ + ri/2-) + C(l + |y|2)e---^ 

If fi'^ is an ergodic invariant probability measure of {Ym{T, .))m, then since <j) is bounded for 
/Lt^-almost any y G H we have by the ergodic Theorem the following convergence when M — )• +00: 

M 



m=l ^ 



M 

To conclude, it remains to choose a initial condition y in this non-empty set, and to use Cesaro 
Lemma on the right-hand side of the estimate. 

We notice that if /i"^ is an invariant probability measure, not necessarily ergodic, having a finite 
moment of order 3, then it is enough to integrate the inequality above with respect to /i'^. 

6. Proof of the estimates 



We need to control the terms given in ()19p . according to the decomposition (jlSp . We recall that 
constants C must be independent from the dimension N and the final time T = mr. 

6.1. Estimate of u{T,x) — E[n(T — r, li)]. The Markov property gives 

^.(r, y) = E[0(y(r, y))] = E[n(r - r, y(r, y))]. 
If < K < 1/2, using Lemma [4. 21 and Proposition 15 . 1 1 we get 

\uiT, y) - E[u{T - T, yi)]| < C{1 + (T - T)-i/2+-)e-^(^--)(E|y(T, y) - Y^W,^_^^)^l\l + 
We can write 

y (r, y)-Yi = (e^^ - Rr)y + [ e(--^)^G(y (s, y))ds - TRrG{y) 

Jo 



+ r e^^-''>^dW{s)-V^RrXl■ 
Jo 
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We use the following properties to estimate the first line in this equality: 

\e'^\c{H) < 1 for s > 0, 

\Rt\c{H) < 1; 

\{-B)-^'^^\\<c\.\- 

therefore the first line in the last expression is almost surely bounded by C{t^/'^~'^ + t)(1 + \y\). 
For the second line, we have 

< CT] 

the last term is controlled in the same way: E| (— i?)~^/^"^'^-v/ri?,-xiP < cr. Therefore we have 

(30) \u{T, x) - E[u{T - T, Yi)] | < C(l + \yf){l + (T - r)-'^/^+^)e-^(T~^)T^/^-\ 

We thus understand that to obtain weak order 1/2 requires to be careful in the estimate. Here 
we used Lemma 14.21 instead of Lemma 15.41 otherwise looking at E|y(T, y) — lip would have not 
been sufficient. The control of the other terms must be done in the same spirit. 



6.2. Estimate of Ofc. We have 



k + l 



ak = K < BY{t) - BrYk, Du{T - t, Y{t)) > dt 



k + l 



E / <{B- Br)Yk, Du{T - t, Y{t)) > dt 



+ E / < B{Y{t) - Yk), Du{T - t, Yit)) > dt 
Jtk 

_ 1 I 2 

6.2.1. Estimate of a\. We use the equality Br — B = tRt-B^. We also decompose a\ using expression 

a}/ = -rE / < RrB^R';y, Du{T - t, Y{t)) > dt 



1,2 



-rE / < RrB'^T R'^-^GiYi),Du{T - t, Y{t)) > dt 

ftk+i '^"i 
rE / < RrB^^/^Y R^^-^xi+i, DuiT - t,Y{t)) > dt; 



1=0 

1,1 , 1,2 , 1,3 
1,1 

The idea is to "share" B^ between different factors - thanks to regularization properties of 
the semi-group {R^)k£fi and to Lemma 14.21 - in order to increase the order of convergence. 
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+ 1, 1 1:1 I 1-2 I 

then aj. = a^' + a^' + 
(1) Estimate of 



Jtk 



thanks to Lemma [3. 2 1 and Proposition 15.11 

By taking expectation, thanks to Lemma 14.31 we have 

m— 1 m— 1 



k=l k=l ■'tk U ^} 



tk " {T-tyi^--' 

^ -dt 



<«i + MV'--/V.(r-<)V.-.- 



and we obtain 

m— 1 

fc = l 



(31) Yl ^ + |y|3)r-(i/2-2«)^V2-2. 



(2) Estimate of a]f' 
First we write 

Using Lemma [3.21 we can prove the fohowing useful inequahty: for r < tq and any A: > 1 

^ 1 1 

^^^^ ^^(ITF^OT^-^''- 

Indeed, 



1 1 „ f^" I 1 



r 



^ (Zr)i- (1 + /xor)''' - ^ i) (1 + /ior)''^ 

poo 1 

/■°° 1 / r 



si V'^log(l + /Uqt) 



Since G is supposed to be bounded, the estimate (I32p yields 

^ ^ (^t)^ (1 + ^0^)"" 

With Lemma 14.31 and Proposition 15.21 we can now write 
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and we get 



m—l 

(33) ^l4'2|<C(l + |y|V/'-"'- 

fc=i 

(3) Estimate of a].'^ 

The analysis of tliis term is more complicated. We recall that since noise is white in space, 
for any t > 0, 'E\{—B)'^W^{t)\'^ < +oo if and only if 7 < 1/4; as a consequence, the strategy 
used above to control a^'^ can not be used directly - otherwise we could only obtain order 
1/4. 

In |;6j, an integration by parts formula is used to deal with the lack of regularity of the 
stochastic integral appearing in the definition of a^''^. An additionnal difficulty arises in 
our situation because the estimate given in Lemma [4.51 is not uniform with respect to time. 
Instead, we remark that the two problems - lack of regularity and bad time dependence - do 
not occur at the same time; therefore a decomposition of the interval [0, tk] into [0, tk — 1] 
and [tk — l,tk] - for k large enough - can help to treat the problems separately. 

Let us explain more concretely the situation at the continuous time level: we have to treat - 
at the finite dimensional approximation level, but with a bound independent from dimension 
- an expression involving B"^ e^''~^^^dW{s). The idea to get rid of this expression is to 
use an integration by parts formula on the whole interval. We can also see that we can do 
this integration by parts only on a subinterval of size independent from t: indeed if t > 1 

ft ft^i pt 

/ e^'-'^'^dWis) = / e'^'-'^^dWis) + / e^'-'^'^ dW{s). 
Jo Jo Jt-i 

The first term is equal to Jq ^ e^*~^~^^^dVF(s), so that thanks to regularization properties 
of the semi-group (e*^) - see Proposition 12.41 - multiplication by B^ is possible - in other 
words we do not require an integration by parts to get a bound independent from the 
dimension; to treat the second term, the lack of regularity remains but can still be treated 
by the integration by parts, with the advantage of involving a smaller interval, where at the 
discrete time level below we can use a uniform control of all quantities, thanks to Lemma 
14.51 Then the same idea of "sharing" B^ can be used again to get order 1/2. 

Let us now explain develop this program for the discrete time situation: by using (jlOp . 
we make the decomposition 

t k—l 

a^'3 = _rE < RrB'^^/T^R^-^xi+i,DuiT - t,Y{t)) > dt 

= -tE < RrB'^R'^-^'dW{s),Du{T -t,Y{t)) > dt 

Jtk Jo 
ftk+i Atk-i)vo 
= -tE < RrB^R';-^'dW{s),DuiT -t,Y{t)) > dt 

Jtk Jo 

- rE / < / ' RrB^R^-^-dW{s),Du{T - t, Y{t)) > dt. 

Jtk J(tfe-l)VO 
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For the first term - which is equal to when < 1 - we use the Cauchy-Schwarz inequahty 
and we directly get 



\E < / RrB^R'^-^'dW{s), Du{T - t, Y{t)) > \ 

Jo 

p(tk-l)VO 

<(E| / RrB'^R^-^'dW{s)\'^)^/\E\Du{T-t,Y{t))\'^)^/'^ 
Jo 



<C{l + \y\')e-^(^-'\ 

thanks to Lemmas 15.41 and 14.31 and to the following inequality - we remark that in the 
integral below tfc-Zs ^ 1^ 



/•(ifc-l)VO /-(tfe-yvO 

E| / RrB'^R^^-^'dW{s)\'^ = \RrB'^R';-^^\l^^H/s 

»{tfe-l)VO 

Tr {RlB^Rl^^-^^^)ds 



Jo 



< 



/•(tfe-l)VO 2 
n TT, \k-l 

/■(tfc-l)VO 1 

-Jo (1 + Mot)'^-'- 

r+oo 1 

< C / —ds 

-Jo (1 + fiory/^ 

<C, 

when r < tq. Then 

rtk+i /■{tfc-i)vo ftk+i 

/ </ RrB^R^~^^dW{s),Du{T -t,Y{t))> dt\<C e-'^(^-*)dt(l + |y|2)r. 

Jtfc Jo Jtk 

For the second term, we use the integration by parts formula of Lemma 14.41 to get 



P'^k + l r^k 

tK < RrB^R^-^=dW{s),Du{T -t,Y{t)) > dt 

Jtk J{tk-l)VO 

r^k+l ftk , s 

= -rE / / Tr(R^-^'B'^RrD'^u{T-t,Y{t))DsY{t))dsdt. 



'tk J(tk-i)yo 

U h e H, {tk - 1) V < s <tk <t < tk+i, with 1^3^ we see that 

D^Y{t) = D^Yk + ['{BrD'^Yk + RrG{Yk)D^Yk)d\ 
+ RrD':{Wit)-W{tk)) 

= D'lYk + {t- tk){BrD^,Yk + R^G{Yk)D'^Yk). 
Therefore, since \tB-j-\c(h) — G 

\D'lY{t)\p<c\D^^Yk\p, 
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and taking supremum over h with \h\ < 1 we get 

(34) \i-BfDsY{t)\ciH) < c\{-BfD,Y,\ciH)- 

The last quantity is estimated thanks to Lemma [4.51 



\D':Y,\^ < C{1 + Lgt)'^'''{1 + \ )\h\. 

When r < To and (tfc — 1) V < s < < t < tfc+i, we see that (1 + Lgt)^^^'' is bounded by 
a constant. 

1 3 

We can then control the second term of a^,' with 
tE f'^' r |i?.(-i?)i/2+2'c|^^^^|(_5)i-3t^fc-;.|^^^^Tr((-B)-i/2-f) 

Jtk J(tfc-l)VO 



X 



f^" .-l+3# 1 1 



X (1 + + (^_')i„J e-^(^"*)(l + 

using Proposition 15.21 and Lemmas 14.51 and 13.21 
On the one hand, we have 

*fc -i+3« 1 /■*''■ 1 1 

ii, ; ^ 77 — TTT^C^S < / , ., ^ , ds < C < +00, 

fe-i)vo '"'^ (l + Mor)('=-'^)^2 "io (1 + /ior)^ " 



for r < To, thanks to (j32 
On the other hand, 



^ /"tfe-l-l 1 -| /'+00 -I 1 

y / (1 + 77^ ^ + 77^ -^)e-^(^^*)dt < / (1 + - + ^)e-^*(it < +00. 



k=l 



Therefore 



(35) Yl < C{1 + \y\')T 



m—l 

2n 1/2-2K 



fc=l 



6.2.2. Estimate of a\. We decompose using the definition oiY - see (I13p : 

af'^ = E / (t-tk) < BBrYk, Du{T - t, Y{t)) > dt 
Jtk 



= E / '^'(t-tfc) < BRrG{Yk),Du{T -t,Y{t)) > dt 
Jtk 

al'^ = E f < I BRrdW{s),Du{T -t,Y{t)) > dt; 



2 2,1 , 2,2 , 2,3 

then al = + a^' + a^' . 
(1) Estimate of a^'"^ 
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Since BB^ = RrB^, a^'^ is bounded by the same expression as a\: by (j3ip . (|33p . (jSSp we 



have 

m— 1 



(36) Yl ^ + + r-(V2-2M)^l/2-2. 

fc=l 

2.2 



(2) Estimate of a^' 

We have 



I 2 



:'\<tE r^'\{-Bfl^+'^Rr\ciH)\G{Yk)\\{-BY/^-'^Du{T-t,Ymdt 
Jtk 



We then have 

m— 1 

(37) I«f I ^ 



fc=i 

(3) Estimate of 

2 3 

We again use the integration by parts formula to rewrite a^' : 

al'^ = E < 1 BR^dW{s),Du{T -t,Y{t)) > dt 



= E j"^^ j TT{RrBD^u{T -t,Y{t))DsY{t))dsdt. 

From (fT3l) . for tk < s < t < t^+i we have D^Y{t) = Rrh; as a consequence, we do not need 
to use the same trick as in the control of a].''^. 
Then we have 



I < E / {t- tk)Tr{RrBD'^u{T - t, Y{t))Rr)dt 
Jtk 

< cr f"^' |i?.(-i?)V2W2|^^^^Tr((-i?)-V2-«/2)|(_5)«^^|^^^^ 

\{-Bf''^^''l'^D'^u{T - t, Y{t)){-Bf'^'^l^\c(^H)dt 

< c(l + |y|2)ri/2-3«/2 1_ \ )e-^(^-^)dt. 



Therefore 

m—l 



(38) ^|af |<C(l + |y|V/'-"^/'- 

fc=i 

With the previous estimates on a}^ and a'^, we get 

m—l 

(39) ^ lofcl < C(l + |y|3)(l + T-(i/2-2«))^i/2 



k=l 
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6.3. Estimate of ftfc- We have 

bk = E ['^' < G{Y{t)) - RrG{Yk),Du{T - t, Y{t)) > dt 

= E / < (/ - Rr)G{Yk), Du{T - t, Y{t)) > dt 

+ E / < G{Y{t)) - G{Yk), Du{T - t, Y{t)) > dt 
Jtk 

■■=bi + bl. 

6.3.1. Estimate ofh\. This term is easy to treat: we have 

|6i|<E r^'|(-i?)-V2+'^(/-i?oi^(H)|G(n)||(-i3)i/2-'^i)n(r-t,y(t))|dt 
Jtk 

<Cri/2- + - l-—)e-i'iT-')dt, 

Jtk {T-ty/^-' 

where we have used Proposition 15. H and the following inequality for < /3 < 1: 

(40) \{-B)-^{I-Rr)\ciH)<Gf,T^ 
Then we see that 

1 

(41) ^ 

k=l 

6.3.2. Estimate of bf.. To estimate we write the scalar product in coordinates with respect to 
the orthonormal basis (fi), and then we expand the terms thanks to the Ito formula. 

If we note Gi =< G, fi > and di =< D., fi >, we have 

< G(y(i)) - GiYk),Du{T - t, Yit)) >= Y,iG^iY{t)) - GiiYkMu{T - t, Yit)). 

i 

The above sum is finite, because we work with finite dimensional approximations. 
Ito formula gives for t^ < t < t^^i 

Gi{Y{t)) - Gi{Yk) = i f Ti{RrRlD^G^{Y{s)))ds 
^ Jtk 

+ I < BrYk,DGi{Y{s)) > ds 
Jtk 

+ I < RrG{Yk),DGi{Y{s)) > ds 
Jtk 

+ [ < DGiiY{s)),RrdWis) > . 
Jtk 

We naturally define bfj'' , for j G {1,2,3,4}, and we now control each term. 

(1) Estimate of b^.'^ 

By definition, we have 

bl^= T'^'eI fy^Tr{R^R;D^Gi(Y{s)))dsdMT-t,Y{t))dt. 
Jtk ^ Jtk i 
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Using the orthonormal basis {fk)k given by assumption 12.21 and recalling that the sums are 
finite, we can calculate: 

^ Tr{RrR*D^G^{Y{s)))dXT - t, Y{t)) = J2 Tr{D^Gi{Y{s))RrK)d^u{T - t, Y{t)) 

i i 

= EE < ^^'^^(n^)) ^^)^. ^2 /^'/^- > d^u{T-t,Y{t)) 
= E E Jy^J^I^^G,{Y{s)).{f,J,)d,u{T - t, Y{t)). 

Using the Cauchy-Schwarz inequality (where j is fixed), we get 
Y,D^G^iY{s)).{fJJJ)^MT-t,Yit))\ 



The second factor of this expression is \(—B)^Du{T — t,Y{t))\H', we control it thanks to 
Proposition 15.11 The first factor is controlled thanks to Assumption 12.51 



since {fj)j is an orthonormal system. 
Therefore 



<C\fj\H\fj\H<C, 



E TriRrR;D^GiiYis)))dMT - t,Ym 

i 

1 °^ 1 

< C{1 + |y|2)(l + ^_)e-'^(^-*) y 



< c(l + + - le-'^^^-^V-i/^-^ (//jr)i/^+" 1 



Then 



and 

m—l 

(42) EI^'''|^^(1 + NV/'-^ 

(2) Estimate of 6^ 



fc=i 

l2,2 
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6f = E 



Thanks to ([9]) and (jlOp . we have 



/ / J2<^-^ry + BrrY,Rr''G{Yi),DGi{Y{s))>diu{T-t,Y{t))dsdt 

Jtk Jtk i l = Q 

+ e/ y2<^^ R^^-^^dW{r),DGi{Y{s))> diu{T -t,Y{t))dsdt 



,2,2,1 , ,2,2,2 



(i) For the first term, recalhng that = BR^ and that G is bounded, we have 
|62.2,i|^|]E I < DG(Y{s)).{BrR';y + BrTY,R';-^G{Yi)),Du{T - t,Y{t)) > dsdt\ 

•'tfe ■Jtk l — Q 

<E / / \{-BrR,\c^H^{\{-B)'-^R'^y\+TY,\i-B)'-^R'r-UH)\m)\) 



'tk Jtk 

X \Du{T -t,Y{t))\dsdt 



< 



if T < To - see (f32]l . 
Therefore 



m— 1 m— 1 ^ „^ 

Jo t 



1 nk+1 



< GT^-^il + \yf)T-'-^/^-^^\ 

(a) For the second term, we again use an integration by parts, after a decomposition 

1 3 

the time interval - as in the estimates for . First, 

52'2.2 I V < 5^ /" ' R''-^-dW{r),DG^{Y{s)) > diu{T - t,Y{t))dsdt 

Jtk Jtk j Jo 

ftk+i ft ^ /■(tfc-l)VO 

= E /zZ<^^/ Rr'^'-dWir), DGi{Y{s)) > diu{T - t, Y{t))dsdt 

Jtk Jtk i Jo 

+ E f'^' f y,<Br f" Rr'^' frn, fj > dp^{r)d,Gi{Y {s))diu{T - t, Y {t))dsdt 
Jtu Jti. ■ ■ _ J(ti.-i)yo 



'tk Jtk ij-^ J{tk-i)yo 

,2,2,2,1 , ,2,2.2,2 

=: ' ' +61.' '. 
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2 2 2 1 

For bfj ' ' , we can work directly and see that 
*fe+i ft ^ /■(tfc-i)vo 



1^2,2,2,1 1 < |]g / / V < 5r / i?^"''-(i^(r),i?G,(y(s)) > ain(r-t,y(t))ds(it| 

Jtfc Jtfc . Jo 

ftk+i ft /•(tfc-i)vo 

< / / E\ < DG(Y{s)).B^ R';-^-dW{r),Du(T-t,Y{t))>\dsdt 

ftk+i ft /-{tfe-yvo 

</ / (E|5r / i?^-''-dTy(r)|2)l/2(]E|^^(2^_i^y(i)) > |2)l/2^^^^ 



<Cr r'"e-A{T-t)(i + |y|2)^ 
thanks to Lemmas 14.31 15-41 and to the following estimate for r < tq 

/■(tfc-l)VO ^(tfc-l)VO 

E\Br / i?^~''-(iVF(r)p < E|b2^^ / R^-^- dW {r)\'^ < C, 

Jo Jo 

1 3 

thanks to the estimate proved to control a^' . 

2 2 2 2 

For bj^' ' ' , we can write thanks to a Malliavin integration by parts and with the chain 
rule 

rtk+i rt . . rtk 



^2,2,2,2 f'^' f y^<Br r > dp^{r)d,Gi{Y{s))diu{T - t,Y{t))dsdt 

J^k Jtk ij^rn ^(tfc-l)VO 

= E f r V < BrR';-'^fm, fj > a|„Gi(y(s)) < Dl'^Yis), fn > d^T - t, Y{t))drdsdt 

Jtk JtkJ{tk-l)^Oij„^^„ 
ftk+1 ft ftk 

+ E / / / V < BrR';-'^ fm, fj > d,Gi{Y{s))df,,u{T - t, Y{t)) < D'^Yit), /„ > drdsdt 

Jtk JtkJ{tk^l)yOi^j^rn,n 

= E f'^' f f" V D^G^{Y{s)){BrR';-'^fm, D^Y{s))diu{T - t, Y{t))drdsdt 
Jtk Jtk J{tk-i)\/o 

+ E / / / ' V < Bi{s, t)BrR'^-^^fm, > drdsdt 

Jtk Jtk J{tk-i)yo 

= ^ / / / Y,^^\iDrYis))*D^Gi{Y{s))BrR'^'^')diu{T-t,Y{t))drdsdt 

Jtk Jtk J{tk-i)yo • ^ ^ 

+ e['^' [ r y2TT({DrY{tyBi{s,t)BrR';~^Adrdsdt, 
Jtu Jti. J(tt~i)yo ■ ^ ' 



'tk Jtk J(tfe-l)VO ' 

where we define a linear operator on H by 

+ 00 

< Biis, t)h, k>=< DGi{Y{s)), h>Y, dln<T - t, fit)) <k,U> 

n=0 

=< DGi{Y{s)), h X D^uiT - t, Yit)).f„ k > . 
We have Y.i < 13i{s,t)h,k >= D^u{T - t,Y{t)).{DG{Y{s)).h, k), and 
\^B,is,t)\ciH) < \DGiYis))\ciH)\D^uiT-t,Yit))\ciHy, 

i 
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so we can write, for (tfc — 1) V < r < 
|^Tr((Z),y(t)*e,(s,t)S,i?^''- 

i 

<\D,.Y{t)\c^H)\Y.^^{s,t)\ciH)\{-BY-'^^^ 

i 

< -1/2-2k.-1+3k/2 1 p-KT-t) 

- (1 + ^or)(^-^-)3«/2 

using Proposition 15.21 Lemma [4. 51 - since {1 + Lqt)^^^^ < C - Lemma [3.21 and estimate (j34p . 
The other term is a little more complicated, because we are not able to control D'^G{Y{s)) 

2 1 

in H. We proceed as in the estimate of 6^' , and we directly calculate the trace. 



I Tr ({DrY{s)YD^Gi{Y{s))BrR^,-'^) diu{T - t, Y{t))\ 

i 

< \DrY{s)\ciH)\Yl Tr [D^Gi{Y{s))BrR';-'^) d,u{T - t, Y{t))\ 

i 

< \DrY{s)\c^H)\{-BfDu{T - t,YmH Y,\{-Br^D^G{Y{s)).U„ f,)\- 



. '(l + ^j-r) 

thanks to the Cauchy-Schwarz inequality. 

By using th 
is bounded by 



2 1 

By using the same analysis as in the estimation of bj , we see that the above expression 



C\DrY{s)\c^H)\i-BrDu{T - t,YmH 



l + k-lr ' 



but the last sum is equal to Tr {Bt-R^ so that we see that indeed the two expressions in 
6^'^'^ are bounded by the same expression. 
Therefore 

IJ.2,2,2,2| 

ftfc+i ft ft)^ —nlT—t) 1 

< C{1 + \yfy I''" {1 + jfhfi^''"'^'"'"- 

as already proved - see (I32p . 

For t.^ 

m— 1 



2221 2222 

Now gathering estimates for bj ' ' and 6, ' ' ' , we obtain 



(43) <^^(i+iyi'y/'"'" 

k=l 
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(3) Estimate of 6^'^ We have 



bl'^ = E 



I Y.< RrG{Yk),DG^{Y{s)) > diu{T-t,Y{t))dsdt 

Jtu J tk 



'tk ■Jtk 



■k + 1 



ftk 

= E I < Du{T -t,Y{t)),DG{Y{s)).{R^G{Yk)) >dsdt. 

Jtk Jtk 

Using that G and DG are bounded, we easily see that 

\bl''\<Cil + \y\')r f"\-^(^-'Ut, 
Jtk 

and that 

m— 1 

(44) 5]|6f|<C7(l + |y|V. 



fc=i 



(4) Estimate of bl''^ 



We use the integration by parts formula of Proposition 14.41 to get 
ftk+i ft 

= DGi{Y{s)),RrdW{s) > d,u{T -t,Y{t))dt 

Jtk Jtk ■ 

= I TT(^{DsY{t))*D^u{T -t,Y{t))DG{Y{s))Rr^dsdt 

Jtk Jtk 

= E f"^' f Tt (^RrD^u{T -t,Y{t))DG{Y{s))Rr^ dsdt, 

using the identity D^Y{t) = Rrh when tk < s <t < tfc+i, as in the estimate of a^'^. 
Now 



X \DMT-t,Y{t))\ciH)Tr{{-B)-^/^-^)dsdt 

< Gil + \y\'y/'-^ £'"(1 + j;f^y-'^^^''>dt, 

and 

TO — 1 

(45) j2 ^ + ivi'y^'"'- 

k=l 

6.3.3. Estimate ofbk-- conclusion. With (gH), (021), (gSD, (glD and (05]), we get 

m— 1 

(46) < Cri/2-2K_ 

fc=i 

6.4. Estimate of c^. We have, using the symmetry of R^, 

h - ^RrR*r = Rr{I " Rr)* + " ^r)(/ " Rr)* , 
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and 

rtk+i 



1 ftk+l 

Ck = -E Tr((/ - RrRl)D'^u{T - t, Y{t)))dt 

= -E / Ti{{I - Rr)il - Rr)*D^u{T - t,Y{t)))dt 
+ E / Tr{Rr{I - Rr)*D'^u{T - t, Y{t)))dt 



6.4.1. Estimate of c\. We have, using inequality (l40l) 
2 At. 

\{-B)y^-^DMT-t,Y{t)){-Bf'^-^\c(H)dt 



•-k 

X 



fc+1 



< C(l + |y|^) / - R,)\c^H)\I - Rr\ciH)'^r{{-B)- 

Jtk 

X (1 + —1— + ^_^)e-A(r-*)dt 
^ {T-t)'^ {T-tY'^' 

Then 

m— 1 

(47) J^|4|<C(l + |y|V/'-'^ 



fc=i 



6.4.2. Estimate of cl. We have, using inequahty (|30 



|c|| < E / Tr((-S)-i/2+-i2,(/ - 

x|(_i5)l/2-«Z)2^(T-t,y(t))(-i3)l/2-.|^^^^rfi 

Jtk 

X (1 + _1_ + \^ )e-^(^-*)dt 

Then 

m— 1 

(48) ^|ci|<C(l + |y|V/'-'^ 

k=l 

6.4.3. Estimate of c^: conclusion. With (|47p and (j48p . we get 

m—l 

(49) ^|cfc|<C7(l + |y|2)rV2-3^ 

k=l 
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6.5. Conclusion. We put together estimates (f39|) . ([16]) . (jl9]) and ([30|) : then passing to the limit 
with respect to dimension, we get the result. 
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